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Abstract 

A recently proposed method for the characterization and analysis of local 
equilibrium states in relativistic quantum field theory is applied to a simple 
model. Within this model states are identified which are locally (but not 
globally) in thermal equilibrium and it is shown that their local thermal 
properties evolve according to macroscopic equations. The largest space- 
time regions in which local equilibrium states can exist are timelike cones. 
Thus, although the model does not describe dissipative effects, such states fix 
in a natural manner a time direction. Moreover, generically they determine 
a distinguished space-time point where a singularity in the temperature 
(a hot bang) must have occurred if local equilibrium prevailed thereafter. 
The results illustrate how the breaking of the time reflection symmetry at 
macroscopic scales manifests itself in a microscopic setting. 



1 Introduction 



Within the framework of relativistic quantum field theory, a general method has 
been established in for the characterization and analysis of states which are 
locally close to thermal equilibrium. This novel approach is based on the idea of 
comparing states locally with the members of a family of thermal reference states 
consisting of mixtures of global equilibrium states: If a given state happens to 
coincide at a space-time point x with some such reference state to some degree 
of accuracy, i.e. if the expectation values of a sufficiently large number of distin- 
guished local observables at x coincide in the two states, the given state is said 
to be locally in equilibrium at x. One can then determine the thermal proper- 
ties of the state at x by noticing that the local observables have an unambiguous 
macroscopic interpretation in all global equilibrium states. Moreover, the space- 
time evolution of these thermal properties is directly linked to the microscopic 
dynamics. 
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The application of this method to concrete models thus requires the following 
steps: (1) Determine the global equilibrium states (KMS states) of the theory; 
their mixtures define the convex set of thermal reference states. (2) Select for each 
space-time point x a set of suitable local observables which are sensitive to thermal 
properties. As was explained in pQ, a natural choice consists of the basic observable 
point fields of the theory as well as of density like quantities which are obtained by 
taking normal products. (3) Determine the thermal interpretation of these local 
observables. This is accomplished by computing their expectation values in the 
global equilibrium states. The result is a family of thermal functions depending 
on the parameters characterizing the thermal states, such as temperature and 
chemical potential; they describe certain specific intensive macroscopic properties 
of the global equilibrium states. (4) After these preparations, one can characterize 
the local equilibrium states in a given space-time region as follows: at each point 
of the region, the mean values of the distinguished local observables taken in such 
a state coincide with those in some thermal reference state. This reference state 
may vary from point to point. One then ascribes to the local equilibrium state 
those macroscopic thermal properties of the respective reference states which are 
fixed by the distinguished thermal functions. The space-time evolution of these 
properties can be derived directly from the underlying microscopic dynamics. 

In the present investigation we apply this scheme to the theory of a free mass- 
less scalar field. In view of the simplicity of this model, it may be necessary to 
comment on the physical significance of such an analysis. It is clear from the outset 
that a non-equilibrium state in this model will never approach (local) equilibrium 
because of the lack of interaction. So we will have nothing to say about this intrigu- 
ing problem. But, appealing to physics, one knows that, once local equilibrium 
has been reached in a system, there appear in general long-living macroscopic 
structures, such as hydrodynamical flows, where dissipative effects are of minor 
importance. Moreover, if the local temperature of a system is sufficiently high and 
the interaction is asymptotically free, these dissipative effects ought to be small 
also on dynamical grounds. Thus one may expect that massless free field theory 
provides some physically meaningful qualitative information about the possible 
space-time patterns of local equilibrium states in hot relativistic systems. 

The present investigation is primarily intended as an illustration of the general 
ideas expounded in P, but it offers some surprises. First, we will see that in 
all local equilibrium states of the model one can determine locally macroscopic 
properties, such as the particle density (distribution function), without having to 
rely on an ad hoc procedure of coarse graining. Second, we will establish evo- 
lution equations for these macroscopic observables, akin to transport equations. 
These equations, together with further constraints arising from the statistical in- 
terpretation of the microscopic theory, imply that the largest space-time regions 
in which local equilibrium states (not being in global equilibrium) can exist have 
the form of timelike cones. Phrased differently, the time reflection symmetry is 
broken in these states, so one can fix with their help a time direction. Finally, it 
turns out that such states determine generically a unique space-time point where 
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a singularity in the temperature (a hot bang) must have occurred if the state was 
in local equilibrium thereafter. An example of such a state was already presented 
in p. The present analysis shows that the appearance of such singularities is 
a generic feature of the model. Incidentally, it provides a consistent qualitative 
picture of the space-time patterns of the cosmic background radiation. 

The general mathematical framework and scheme for the identification and 
analysis of local equilibrium states is outlined in the subsequent section. For a more 
detailed discussion and physical justification of the method we refer to In Sec. 3 
we describe the model and carry out the preparatory steps, indicated above, which 
are required for the analysis of the local equilibrium states. Section 4 contains the 
discussion of the space-time behaviour of thermal properties of these states and 
the derivation of macroscopic equations governing their evolution. In Sec. 5 we 
show that local (in contrast to global) equilibrium cannot exist in all of space- 
time and we determine the maximal regions admitting this feature. The paper 
concludes with a discussion of these results and some interesting perspectives. 



2 General framework and scheme 

The basic objects in our analysis are local observable fields, such as conserved 
currents and the stress energy tensor, denoted generically by (f>(x) (with tensor 
indices omitted). Mathematically, these fields may be interpreted in different ways. 
Usually, they are regarded as operator valued distributions, i.e. their space-time 
averages with test functions / having compact support in Minkowski space M 4 , 

0(/)= fdxf(x)<j>(x), (2.1) 

are assumed to be defined on some common dense and stable domain T> in the 
underlying Hilbert space Ti. The finite sums and products of the averaged field 
operators generate a *-algebra A, called the algebra of observables, consisting of 
all polynomials A of the form 

A = J2<P(fMh) ■■■<!>( fn)- (2-2) 

We include also all multiples of the unit operator 1 in A. As the fields are real 
(being observable), the *-operation (Hermitian conjugation) on A is defined by 

where / denotes the complex conjugate of /. The Poincare transformations A G V\ 
act on A by automorphisms a a, 

a x (A) = 0(/i,a)0(/ 2 ,a) • ■ • <P(fn,x), (2-4) 

where f\ is defined by f\(x) = D(X)f(X^ 1 x) and D is a matrix representation of 
V]_ corresponding to the tensor character of the respective field <p. If A is a pure 
translation, A = (1, a), a G R 4 , we denote the corresponding automorphism by a a . 
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The states of the physical systems are described by positive normalized expecta- 
tion functionals on the algebra A, generically denoted by u. We recall that these 
functional have the defining properties oj{c\Ai + c 2 A 2 ) = Ciu(Ai) + c 2 cg>(A 2 ), 
Ci,c 2 G C (linearity), u(A*A) > (positivity), and lu(1) = 1 (normalization). As 
has been explained in the introduction, an important ingredient in the present in- 
vestigation are the global equilibrium states which are characterized by the KMS 
condition [2 Ej. Given a Lorentz frame, fixed by some positive time-like vector 
e G V + of unit length, e 2 = 1, this condition can be stated as follows. 

KMS condition: A state up satisfies the KMS condition at inverse temperature 
(3 > in the given Lorentz frame if for each pair of operators A, A' G A there is 
some function h which is analytic in the strip Sp = {z G C : < Imz < j3} and 
continuous at the boundaries such that 

h{t)=up{A'a te {A)), h{t + i/3)=u;p{a te {A)A'), teR. (2.5) 

In this situation ujp is called a KMS state. 

Any KMS state uop describes an ensemble which is in thermal equilibrium in the 
distinguished Lorentz frame, describing the rest system of the state. As we have 
to keep track of both, temperatures and rest systems, we combine this information 
into a four vector f3e G V+, again denoted by f3. Moreover, in order to simplify 
the subsequent discussion, we assume that for any given j3 G V + the corresponding 
KMS state up is unique. As a consequence of this uniqueness assumption, the 
KMS states transform under Poincare transformations according to 

cup o o^ 1 = UJ A/ 3, (2.6) 

i.e. they are isotropic in their rest systems and invariant under space-time transla- 
tions. Next, let us introduce the thermal reference states. They are defined, for any 
compact subset B C V + , as arbitrary mixtures of KMS states cup with (3 G B. We 
emphasize that the restriction to mixtures with compact temperature support B 
is merely a matter of mathematical convenience in the present context. In general 
one may admit larger families of reference states with non-compact temperature 
support. In order to obtain a more explicit description of the reference states and 
to avoid mathematical subtleties, we assume that the KMS states ujp are weakly 
continuous in /3, i.e. all functions 

(3 i-> up(A), A e A, (2.7) 

are continuous. Apart from phase transition points (which are excluded here by 
our uniqueness assumption), this property is expected to hold quite generally. 
With this input, the states ojb G Cb can be represented in the form 

u B {A) = J dp((3) ujp(A), A e A, (2.8) 

where p is a positive normalized measure which has support in B. The convex set 
of all these states will be denoted by Cb and the elements of C = [J Cb, where the 
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union extends over all compact subsets B C V + , will be our reference states for 
the characterization and analysis of the local equilibrium states. 

It is mathematically meaningful and conceptually convenient in the present 
analysis to consider also the unregularized fields <p(x). Because of their singular 
nature they have to be interpreted in the form sense, however. Let us illustrate 
this for the KMS states. In view of the invariance of these states under space-time 
translations, the expectation values (</>(/)) do not depend on the choice of / as 
long as the value of the space-time integral of the test function is kept fixed. One 
may thus extend the thermal reference states to the unsmeared fields <f>(x) by the 
formula 

Mfa)) = M<Kf)), ( 2 - 9 ) 

where the space-time integral of / is put equal to 1. This expression does not 
depend on x. We denote the linear space generated by all fields <f>(x) at x by Q x 
and assume without further mentioning that all states to of interest here can be 
extended to the spaces Q x for x varying in some region. 

It is crucial in the present approach that the local fields <p(x) provide the same 
information about properties of the thermal reference states as certain macroscopic 
(central) observables. To verify this, we pick a test function / of compact support 
which integrates to 1 and put f n (x) = n~ 4 /(n _1 a; — x n ), n G N, where x n is a 
sequence of translations tending sufficiently rapidly to spacelike infinity such that 
the support of f n lies in the causal complement of any given bounded region for 
almost all n. It follows from the uniqueness of the KMS states and their invariance 
under space-time translations by an application of the mean ergodic theorem that 
in the limit of large n 

M<t>(f*) Wn)) - M<KUY) M<P(fn)) - o. (2.10) 

Moreover, because of the local commutativity of observables and the support prop- 
erties of the test functions, 4>(f n ) commutes for almost all n with any given A 6 A, 
i.e. the operators <p(f n ) form a central sequence in A. Applying standard argu- 
ments, it follows that the limit 

$ = lim 0(/„) (2.11) 

exists in all thermal reference states and defines a macroscopic observable com- 
muting with all elements A e A and the translations. In fact, making use of (|2.8jl . 
(|2.1()|) and applying the dominated convergence theorem, one obtains 

uj b (A*$A) = lim u B {A*(j)(f n )A) = [ dp(P)up(A*A)u0((l>(O)). (2.12) 

This relation implies that the (according to (|2.7j) continuous) function 

P» $(/?)= ^(0(0)), (2.13) 

called thermal function, determines the central decomposition of $; hence we may 
identify $ with this function. Such central operators are called macro-observables 
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in the following. Since ujb{4>{x)) = u B ((p(0)) = uj b (&), the expectation values of 
the local observables 4>{x) in the thermal reference states can be interpreted in 
terms of the macro-observables $, as claimed. 

Let us now turn to the characterisation of the local equilibrium states. Fixing 
some suitable subspace S x C Q x , we say a state uj is ^-compatible with a thermal 
interpretation at x (^-thermal, for short) if there exists some ujb G C such that 

uj(<f>(x)) = u B ((f>(x)), 4>{x) G S x . (2.14) 

Thus uj cannot be distinguished from the thermal reference state ujb by the local 
observables in S x . One can therefore consistently assign to it the thermal properties 
of ujb described by the macro-observables $ corresponding to <p(x) G S x . More 
explicitly, taking into account that uj((f>(x)) = ujb(4>(x)) = one can lift uj 

at x to the space of macro-observables, setting 

uj{$){x) = cu((f>(x)), <f>{x) G S x . (2.15) 

In this way the state uj acquires at x a physically meaningful (partial) interpretation 
in macroscopic terms. 

This definition extends in a straightforward manner to states which are locally in 
equilibrium in (open) regions Oct 4 . One first identifies the spaces of observables 
S x , x G R 4 , with the help of the automorphic action of the translations, setting 

S x = a x (S ). (2.16) 

This being understood, a state uj is said to be iS^-thermal if for each x G O 
there is some thermal reference state GC such that relation ()2.14|) holds. The 
functions 

x^uj($)(x) =uj(<p(x)), 0(0) eS , (2.17) 

then describe the space-time behaviour of the mean values of the macro-observables 
$ in the state uj. Hence they provide a link between the microscopic dynamics 
and the evolution of macroscopic thermal properties, i.e. the thermodynamics of 
the state. It is this aspect which will be studied in the framework of our model. 

In order to exclude certain pathologies, we restrict attention here to iSo-thermal 
states uj admitting reference states uj Bx G C, x G O, which are weakly integrable 
in x and have temperature supports B x contained in some compact set B C V+ 
for x varying in any given compact subset of O. So the functions x \— > uj(<&)(x) 
are differentiable in the sense of distributions. As a matter of fact, for any test 
function / having compact support in O 

dxf(x)uj($)(x)\ < WfWjng, (2.18) 

where \\fWi denotes the Z^-norm of /, and we have introduced on the thermal 
functions (|2.13J) the seminorms 

||$|| B = sup|$09)|. (2.19) 

This statement follows from relations (|2.8jh ()2.13|) and (|2.14j) . It will allow us to 
extend the lifts of local equilibrium states to larger spaces of macro-observables. 
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3 The model 



In the remainder of this article we apply the preceding general scheme to the theory 
of a free massless scalar field. After a brief synopsis of the model, we proceed 
according to the steps outlined in the introduction and determine the thermal 
reference states as well as a distinguished space of local observables needed for the 
identification and macroscopic interpretation of the local equilibrium states. An 
example of such a state is given at the end of this section. 

The free massless scalar field 4>q{x) on M 4 is characterized by its field equation 
and commutation relation 

□ s o (x) = 0, [M^i),Mx 2 )} = (2tt)- 3 [ dpe~^-^ p e( Po )S( P 2 ) • 1, (3.1) 



which are to be understood in the sense of distributions. It generates, together 
with its normal products with reference to the vacuum state (cf. the examples 
below), a polynomial *-algebra A. This algebra is stable under the actions of the 
Poincare group V]_, given by ot^ a ((j) (x)) = </>o(Ax + a) and the gauge group Z 2 , 
given by 7O0O)) = ~M X )- 

We restrict attention here to states uo on A which are gauge invariant, i.e. 
= u, so the respective n-point functions of 0o vanish if n is odd. The simplest 
examples of this kind are quasifree states. They are completely determined by their 
two-point functions through the formula 

U;(0oOl)0o(>2) • • -0oOn)) 

T, P ^ inss ^(Mxh)Mxi 2 ))---^(Mxi n - 1 )M x iJ) neven, (3.2) 
n odd. 

It is a well known fact that the algebra A has a unique gauge invariant KMS state 
ujp for each temperature vector (3 G V+. This state is quasifree, so it is determined 
by its two-point function given by 

MMxi)Mx2)) = (27r)- 3 |^e-^-^)%(p )5(p 2 ) T 3^. (3.3) 

The KMS states up comply with our continuity assumption ()2.7|) and fix the 
convex set C of thermal reference states which enters into our definition of local 
equilibrium states. 

Next, we describe the spaces S x of local observables which will be used to 
analyze the local equilibrium states. They are generated by the unit operator 1, 
the field <fio{x) and balanced derivatives of its normal ordered square. Introducing 
the multi-index notation 11 = (u 1; u 2 , . . ■ U m ) and setting df 1 = d t ■ ■ ■ d t , the 
latter observables are given by 



BUx) = lim %{fa{x + ()Mx - - ^oo{Mx + 0M% - 0) 1), (3.4) 

7 



where the limit is approached from spacelike directions ( and denotes the 
vacuum (i.e. (3 / , oo) state. We assume that all states of interest here are in the 
domains of these forms and note that 8^(x) = if m is odd since 0o(^ + C)0o( a; — C) 
is even in ( as a consequence of locality. 

As has been explained in Sec. 2, each 9^(x) determines some macro-observable 
QV. The corresponding thermal functions have been computed in pQ and are given 
by 

P _> 0M(0) = ^(^(0)) = c m dg (P 2 )-\ (3-5) 

where c m = if m is odd and c m = (— l) m//2 (47r) m (m + 2)! _1 S m+ 2 if m is even, 
_B n being the (modulus of the) Bernoulli numbers. As dp (f3 2 )~ l = on V + , the 
observables G^ generate only a subspace of the space of all macro-observables. 
Yet, as we will see, the linear combinations of the corresponding thermal functions 
are dense in the space of smooth solutions of the wave equation on V + with respect 
to the topology induced by the seminorms ()2.19|) . This fact is of relevance since 
it implies that many macro-observables of interest, such as the entropy density 
PP or the particle density, can be approximated in the states considered here with 
arbitrary precision by linear combinations of the 0^. 

Lemma 3.1 Let H be any smooth solution of the wave equation on V + , let B C V+ 
be compact and let e > 0. There are finitely many constants c^, c'„ such that 

||E-5> M e"|| fl <e and II^H-^e^|| B < £ , 

where z//x denotes the multi index (z/, /xi, /x 2 , • • • /i m ), v = 0, 1, 2, 3. 

Proof. It is sufficient to establish the statement for closed double cones B of the 
form B — {(k, 0) + V+} fl 0) — V + } for < k < 1. These regions are stable 

under the involution l((3) = j3//3 2 , (3 6 V+. As is well known, this transformation 
induces an involution on the smooth solutions of the wave equation on V + , given 
by J(H)(/3) = t,(j3) 2 H(t(/3)). It is continuous relative to the seminorms (|2.19|) . 

||J(S)|| B < k- 2 ||5|| B , ||^/(H)|| B < 2k-' ||S|| fl + 3/s;- 4 ^ ||^'S|| B . 

(Similar bounds can be established also for higher derivatives by elementary com- 
putations.) Since / is an involution, H and J(H) may be interchanged in these 
inequalities. 

For the proof of the first half of the statement it is convenient to proceed, 
in a first step, to the functions l^Q^, where 1^ = 1^... l^ m , I being any positive 
lightlike vector, and, in a second step, to their images 1(1^0^) under the involution 
/. According to the preceding remarks it suffices to show that the linear span of the 
latter functions is dense in the space of all smooth solutions of the wave equation 
with regard to the seminorms || • The functions 1(1^ G^), resulting from (|3.5|) . 
have the simple form 

/?^/(^)(/?)=4(z/?r\ 
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where c' m ^ if m is even and c' m = if m is odd. So we have to show that suitable 
linear combinations of the functions (5 \— > (//3) m for even m and lightlike / approx- 
imate any given smooth solution of the wave equation, uniformly on compact sets 
B C V+. 

We make use of the fact that the functions z i— > e ±l ^ are analytic in the 
complex half-plane C + = {z e C : Rez > 0}. Thus they can be represented on 
any compact subset C C C + by a uniformly convergent power series in z, whence 
z h e ±12 = gi^v 7 ^ can k e re p resen ted by a uniformly convergent power series in 
z 2 G C. As k\1\ <lf3< 2/t~ 1 |Z| for (3 G B and / = (|Z|, Z), we see that the functions 
/5 i — y e ±ll l 3 can be represented by uniformly convergent power series involving only 
even powers (l/3) m , j3 G B. Moreover, the functions 

[d 3 lf ± (l)e ±il ?, 

where f± are absolutely integrable, can be approximated by linear combinations of 
(5 I— > e ±li ^, uniformly on compact sets -B. So we conclude that linear combinations 
of /? i — > (Z/3) m for even m and suitable / approximate these functions as well. 

It remains to show that the restriction of any smooth solution S of the wave 
equation to the double cone B admits an integral representation as given above. 
This can be seen if one extends the smooth Cauchy data of S on the base of B to 
test functions with compact support in a slightly larger region. Using these data 
as initial values, one obtains a smooth solution Eb of the wave equation which 
coincides with S on B according to standard uniqueness results. In view of the 
regularity properties of its Cauchy data, H# admits an integral representation as 
given above. In fact, the corresponding functions f± are, apart from an integrable 
singularity at the origin, smooth and rapidly decreasing. This completes the proof 
of the first part of the statement. 

The proof of the second part proceeds along similar lines and we therefore 
indicate only the essential steps. We begin by noting that 

where c' m+1 ^ if m is odd. Applying the involution I to f3 i— > (Z/3) m (/3 2 ) _m_1 , 
we see that it suffices to show that suitable linear combinations of the functions 
j3 i— > (l/3) m for odd m and lightlike / approximate any given smooth solution of the 
wave equation, uniformly on compact sets B C V+, and that the same holds true 
for their respective partial derivatives. 

We now use the fact that the functions z i— > (i/Vz) e ±l ^ are analytic on C+. 
Hence z \— > e ±lz = (z/Vz?) e ±l ^ can be represented, on the domains specified 
above, by uniformly convergent power series involving only odd powers of z. So 
the functions (5 \— > e ±ll/3 can likewise be represented by uniformly convergent 
power series involving only odd powers (//3) m , /3 G B. In view of the analyticity 
of these functions, this uniformity of convergence holds also for their respective 
partial derivatives. 
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Now if the functions f± in the above integral representation of solutions of the 
wave equation are sufficiently rapidly decreasing, these integrals can not only be 
approximated by suitable linear combinations of the functions (3 \— > e ±tl @, uni- 
formly on B, but the same holds true also for their respective partial derivatives. 
Yet, as has been explained, the restrictions of smooth solutions of the wave equa- 
tion to B can be represented in this desired way, so the proof of the statement is 
complete. □ 

In order to illustrate the significance of this result, let us consider the mean 
density of particles of momentum p. From a macroscopic point of view, the cor- 
responding observable consists, for a system in a finite box, of the product of the 
annihilation and creation operators of a particle of momentum p, divided by the 
volume of the box. Being in a setting describing the thermodynamic limit, this 
observable is obtained as limit of sequences of elements of A which are bilinear in 
the basic field 4> , 

N p = lim 2\p\ 0o(/p,n)*0o(/p,n)- ( 3 -6) 

n 

Here p = (\p\,p) and f p , n (%) = n~ 5 ^ 2 g(n~ 1 xo) h(n~ 1 x) e ipx , where g, h are test 
functions satisfying jdxog(xo) = 1 and jd 3 x \h(x)\ 2 = 1. The limit exists in all 
thermal reference states and its expectation values can be centrally decomposed 
similarly to relation (|2.12|) . Thus each N p is a macro-observable commuting with 
all elements of A. The corresponding thermal function is, as expected, given by 
Planck's famous formula for the density of massless particles of momentum p in a 
thermal equilibrium state, 

(3 » N p {/3) = (2tt)- 3 — Lp (3.7) 

As these functions are smooth solutions of the wave equation on V+, we infer 
from the preceding lemma that, instead of relying on (|3.6jh one can determine the 
particle density N p also locally. More precisely, for given temperature support B 
there are local observables in S x which allow to determine the mean values of N p 
in all thermal reference states in Cb with arbitrary precision. Since all powers of 
N p are represented by smooth solutions of the wave equation as well, the same 
holds true also for all moments of N p in these states. 

Having thus determined all relevant quantities for the examination of local 
equilibrium states, let us turn now to their actual analysis. We will consider states 
uj which are iS^-thermal in regions O of Minkowski space, where S x , x 6 O, are 
the spaces of local observables defined above. Taking into account the preceding 
lemma and relation (|2.18|) . these states can be lifted at each x £ O to all macro- 
observables H whose thermal functions are smooth solutions of the wave equation 
on V + . These observables will be called admissible in the following. We are 
interested in the space-time behaviour of their respective expectation values, 

x^u(E)(x), xeO, (3.8) 



which describe the evolution of these states from a macroscopic point of view. 
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An interesting example of a local equilibrium state which is Sy + -thermal was 
presented in |T]. This state is quasifree and therefore fixed by its two-point function 
which, for xi,x% G V+, has the form 

^ hb (0o(*i)<M* 2 )) = (2n)~ 3 Jdpe^'^e(po)5(p 2 ) 1 _ e _^ 1+a!a)p , (3.9) 

where r\ > is some parameter. As was shown in pQ, the restriction of u bhb to 5^ 
coincides with the KMS-state UJ/3( X ), where f3(x) = 2r]x, x G V+. So one obtains 
for the expectation values of the admissible macro-observables S 

w bhb (S)(z) = S(2^), (3.10) 

where /3 i— > H(/3) is the thermal function corresponding to S. Plugging the observ- 
able T 2 determining the square of the temperature into this equation (which is 
legitimate since (3 \— > T 2 {(3) = (/5 2 ) -1 is a smooth solution of the wave equation), 
one finds that the temperature of the state u bhb tends to infinity at the boundary 
of V + . As a matter of fact, the state u bhb describes the space-time evolution of 
a hot bang at the origin of Minkowski space. This can be seen most clearly by 
looking at the expectation values of the density operators N p , 

u> bhh (N P )(x) = (27T)- 3 e J p _ v x e V+. (3.11) 

If x approaches the boundary of the lightcone V + , this expression stays bounded 
unless x and the four-momentum p become parallel. Thus the bulk of the particles 
giving rise to the infinite temperature at the boundary of the lightcone originates 
from its apex. Whence the state u bhb provides a qualitative picture of the space- 
time patterns of the radiation caused by a hot bang. We will see that the occurrence 
of such singularities is a generic feature of the present model. 



4 Evolution equations 

We want to demonstrate now how the microscopic dynamics leads to evolution 
equations for the expectation values of macro-observables in local equilibrium 
states in the present simple model. We begin by noting that the field equation for 
0o entails the following equations (in the sense of distributions) 

d xv d <: v Mz + OMx-0 = o, (4-1) 

□* 0o(^ + C)0o(^ - C) = - □< 0o(x + C)0o(x - 0- (4-2) 
This implies for the local observables 8^, introduced in relation (|3.4|) . 

d l/ 6 v t i {x) = 0, (4.3) 

ue^{x) = -e v ^{x). (4.4) 
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In particular, each 9^ is a conserved tensorial current. The latter relations lead to 
linear evolution equations for the expectation values of the corresponding macro- 
observables in local equilibrium states. As a matter of fact, the following 
lemma holds. 

Lemma 4.1 Let uj be an So-thermal state and let H be any admissible macro- 
observable. Then 

d v uj(d v E){x) = 0, □w(S)(z)=0 forxeO, 

where d v E is the macro-observable with thermal function (3 i— > d v, E{(3). 

Proof. For the proof of the first equation we proceed from relation ()4.3j) which 
implies for any choice of constants 

As uj is iSo-thermal, we also have for x £ O 

^Q2^^)(x) = u;Q2^e^(x)), 

and combining these two equations we arrive at 

d v uQ2c lt e"i t ){x) = o. 

The first part of the statement now follows from the continuity property 1)2.18)1 of 
local equilibrium states and Lemma f3. 11 

For the second equality we make use of relation (J4.4j) which implies 

(3 h- QT((3) = u)p{9^{x)) = -u p {U6*{x)) = -Uu p {9^{x)) = 0, 

where the latter equality follows from the invariance of the KMS states under 
space-time translations. Hence Qp = 0, so taking into account once more that uj 
is iSo-thermal we obtain 

= - W (53c^"( a :)) = -oj(J2^T)(x) = 0. 

As before, the second assertion then follows from the continuity properties of uj 
and Lemma fa. II □ 

A macro-observable of particular physical interest is the particle density N p . 
Its evolution in local equilibrium states is described in the following proposition. 

Proposition 4.2 Let uj be any So-thermal state and let p be any positive lightlike 
vector. Then 

(pd)uj(N p )(x) = 0, Duj(N p )(x) = forxeO. 
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Proof. Consider the function L p on V+, 

/3^L p (/5) = (27r)- 3 ln(l-e^). 

It is a smooth solution of the wave equation and consequently L p is an admissible 
macro-observable. Now 

d v L p {fi) = V V (2tt)- 3 -^—^ = p»NM- 
Applying the preceding lemma, we conclude that 

(pd)u(N p )(x) = d v u{p v N p ){x) = d v u{d v L p ){x) = 0, 
as claimed. The second part of the statement follows directly from the lemma. □ 

The preceding result shows that the local equilibrium situations, described by 
iS^-thermal states, are compatible with those considered in the context of trans- 
port equations, such as the Boltzmann equation for the particle distribution func- 
tion h4 N(x,p). In the latter case one characterizes the distribution functions 
N describing local equilibrium by the condition that the collision term in the 
Boltzmann equation vanishes. In the present approach these distributions are to 
be identified with i,p h u(N p )(x), and the first part of the proposition shows 
that these functions satisfy the collisionless Boltzmann equation as well. From 
the second part we see that they have to comply with an additional constraint of 
dynamical origin. 

We mention as an aside that the functions i,p h N(x,p) contain all macro- 
scopic information about the local equilibrium states. As a matter of fact, given 
any non-negative function N satisfying in some convex space-time region O the 
two equations in the proposition and being sufficiently well-behaved in p, one can 
define a quasifree functional (p on the algebra A of local observables, setting for 
xi,x 2 e O 

(p (0o M 00 (X 2 )) 

= ^oo(0o(^i)0o(^2)) + J^pIp]' 1 N((x 1 +x 2 )/2,p) cos((xi -x 2 )p). (4.5) 

This functional is real and satisfies all linear constraints imposed by the field O , 
but it need not be positive. Yet if it complies with the latter condition, it defines a 
local equilibrium state in the region O. As we shall see in the subsequent section, 
the condition of positivity imposes further stringent constraints. 

5 Regions of local equilibrium 

In this section we want to determine the shape of space-time regions in which 
local equilibrium is possible. We restrict attention here to convex regions. Given 
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a state u, let O w be any maximal (in-extendible) region of this kind in which uj 
is Sq u — thermal. As was exemplified in Sec. 3, there are states where O w contains 
a lightcone. It will turn out that, conversely, any such O u has to be contained 
in some timelike cone unless the lifts of u to the admissible macro-observables 
coincide at all points of O u with some fixed thermal reference state. 

It will be convenient in this analysis to proceed from uj to suitably regularized 
states. To this end we pick some non-negative test function / which has support 
in a small neighbourhood of £ M 4 and satisfies Jdy f(y) = 1. The corresponding 
regularized state is given by Uf = fdyf(y) u>oa y , hence w/oa I = w/_ I . As the set 
C of thermal reference states is convex, Uf is (So-thermal in some slightly smaller 
region O C O w . Thus, taking into account relation 1)2.18)1 and Lemma 13.1) we 
obtain for all admissible macro-observables S and multi indices /x the bounds 

\d»u f {E){x)\ <c"||S||b, (5.1) 

uniformly on compact subsets of O, where c^ are certain constants and B C V + 
is compact. 

We will analyze the states ujf on subspaces Y p of macro-observables, fixed by 
the positive lightlike vectors p. These spaces consist of operators whose thermal 
functions are of the form V+ 3 [3 i— > T(/3p), where T is smooth on WL + . So they 
are admissible macro-observables. As their sums, products and adjoints are also 
of this form, the spaces T p are in fact abelian *-algebras. Of special interest are 
their elements E p and M p given by 

(3 E p {0) = e^ p , (3 » M p n (f3) = {(3p)\ n e N„. (5.2) 

Since uif is a state, its lift u)f( - )(x) to the admissible macro-observables at any 
point x £ O defines a state on each T p . We collect some pertinent properties of 
these lifts in the subsequent lemma. 

Lemma 5.1 Let ujf be any regularized So -thermal state. Then, for x £ O, 

(a) (pd)u f (E p )(x) = 0, nu f (E p )(x) = 0, 

(b) \uj f (E p )(x)\< 1, 

(c) p i— > u)f(E p )(x) extends to an entire analytic function on C 4 which is smooth 
in x, 

(d) < LOf{M p n+l ){x) < (^(M;)(x)) n/{n+1) {u f (M p ^)(x)) 1/{n+1 \ n £ N . 

Proof. Statement (a) follows from Lemma (j4.1j) and the relation p v E p = —i d h ' E p . 
For the proof of (b) we make use of the fact that Uf{ ■ ){x) is a state on Y p , hence 
\ujf(E p )(x) | 2 < u>f(E p *E p )(x) = 1, x £ O. In order to establish the existence of 
the extension in (c), we first extend the operators E p to admissible operators Ek, 
k £ C 4 . Their respective thermal functions are given in proper coordinates by 

(3^E k {{3) = (cos(/3 |fc|)+«A; |fc|- 1 sin(/3o|fc|))e-^. 
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These functions are entire analytic in k and locally uniformly bounded if (3 varies in 
compact sets B C V + . It thus follows from the continuity properties subsumed in 
(|5.1|) that k i — > ujf(Ek)(x) provides the desired extension. The remaining statement 
(d) is a straightforward consequence of Holder's inequality, taking into account 
that the operators M p n , n E N , are positive elements of T p (being represented by 
positive functions on V+). □ 

After these preparations we can establish now the limitations on the possible 
shape of regions of local equilibrium, indicated above. 

Proposition 5.2 Let uj be a state which is So^-thermal in some convex region 
containing a lightcone. There are the following alternatives for the lifts uj{ ■ )(x) of 
uj to the admissible macro-observables: 

(a) The lifts do not depend on x £ O^, i.e. they all coincide with some fixed 
thermal reference state. 

(b) The lifts depend non-trivially on x £ O u , and there is some timelike simpli- 
cial cone (an intersection of characteristic half spaces) containing O u . 



Proof. As explained above, we proceed from uo to the regularized states Uf. Any 
one of the latter states is iS^-thermal in the region O = f] x g s /{Cw + %} , which 
likewise is convex and contains some lightcone, say V + for concreteness. Because 
of the convexity of O, V + + y C O for y £ O. We consider now the functions, p 
being an arbitrary positive lightlike vector, 

x I— > E p {x) = u)f(E p )(x), x £ O. 

Introducing the coordinates x± — Xq ± ex and x± = x — (ex)e, where e = p/\p\, 
it follows from part (a) of the preceding lemma that E p does not depend on x + 
and satisfies Aj_E p (x) = 0, x £ O, where Aj_ is the Laplacian with respect to x±. 

Now, given x_ £ {y_ : y £ O} and any x± £ M 2 , there exists an x £ O 
with these components. For if y £ O, the point x with the components x_ = y_, 
x + = y + + t, x± is, for sufficiently large t, contained in the lightcone z + V + C O, 
provided z £ O and (y — z) is positive timelike. So for fixed X- £ {y_ : y £ O} 
one has A±E p (x) = 0, x± £ M 2 . But according to parts (b) and (c) of the 
preceding lemma, x ^ E p (x) is smooth and bounded in modulus by 1, hence it 
cannot depend on x± in view of the growth properties of non-trivial solutions of 
the Laplace equation (Harnack's inequality). 

Because of the analyticity properties of E p established in part (c) of the lemma, 
we can represent this function as power series, E p (x) = J2m=o \p\™ c ^i x ) e/X ' wnere 
e = (l,e) and all coefficients c^(x) are smooth in x. Moreover, in view of the 
results obtained in the preceding step, x \— > c^{x)e^ can depend only on x__ = ex 
in a non-trivial manner. Thus, by /c-fold partial differentiation with respect to x 
and an application of the chain rule, we obtain 

M)%(2)e^ = (ye) k d k CtJl {x)e», x £ O, y £ R 4 . 
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Expressing the products of the components of e in terms of spherical harmonics, it 
is apparent that, for k > m, this equality can only be satisfied if <9 fc c^,{x)e^ = 0. 
So each function x \— > c^(x)e^ is a polynomial in ex of degree m or less with 
coefficients which are polynomials in the components of e. 

Making use of the fact that E p (x) = J2Z=o ^Sr" w /( M e w )0), where M e m are 
the operators introduced in (|5.2|) . we therefore find that each function 

,x i-> W/ (M e m )(x) = (-z) m m! c M (x) e^, x G 0, 

is a polynomial in ex of degree m or less. According to part (d) of the preceding 
lemma these polynomials are non-negative. Moreover, there is some m such that 
the corresponding polynomial is of first degree in ex, unless x t— > ojf(M™)(x) 
is constant for all m G N. To verify this, let x i— > cj/(M e fc )(x) be constant for 
k = 0, . . . m — 1. It then follows from (d) and the fact that the polynomial degree 
of x h- > a;/(M e 2m )(x) is at most 2m — 1 that x \— > cj/(M e m )(x) increases no 
faster than l^ 2 ™- 1 )/™- f or large x G V+ C 0. Hence x i— > u;/(M e m )(x), being a 
polynomial, has to be of first degree in ex if it is not constant. 

If the second alternative in the statement holds, i.e. if x \— > cj/(S)(x) depends 
non-trivially on x for some S, there is also some M e m for which this is true since 
the linear span of the latter operators is dense in the space of admissible macro- 
observables, cf. the proof of Lemma 13.11 Thus, because of the preceding results, 
we may assume that 

cu f (M e m )(x) =P f (e)(ex)+Q f (e), x G O, 

where P/, Qf are polynomials and Pj is not identically zero. But u;/(M e m )(x) > 0, 
x G 0, so 

P f (e) (ex) +Q f {e) > 0, x G 0. 

Since V+ C and Py 7^ it follows that P/(e) > for almost all e and Qf(e) > 0. 
Hence C f] e {x : ex > —Qf(e)/Pf(e)}. In particular, is contained in some 
positive timelike simplicial cone consisting of the intersection of four characteristic 
half-spaces. The same is therefore true for O w since the support of / is compact. 

We mention as an aside that, due to the linearity of Uf in / and Uf oa y = uif_ v , 
one has P f (e) = Jdyf(y)P(e) and Q f (e) = Jdy f(y) {(ey)P(e) + Q(e)) for 
certain polynomials P, Q which do not depend on /. 

The statement has thus been established for regions 0^ containing V + . But the 
arguments can be carried over to arbitrary regions containing some shifted forward 
or backward lightcone, so the proof of the proposition is complete. □ 

The preceding result shows that non-trivial local equilibrium states u fix a 
time direction in the sense that they determine some maximal timelike cone in 
which local equilibrium can exist. In spite of the fact that the underlying theory 
is time reflection invariant, an observer who has such a state as a background may 
consistently take the corresponding cone as an indicator of the future direction. 



16 



If he perturbs uo by some local unitary operation U, u>jj( ■ ) = uj(U* ■ U), he finds 
that the perturbed state ujj coincides with u with regard to all measurements 
performed after the perturbation, i.e. in the future of the space-time support of 
U. This is a consequence of the timelike commutativity of observables (Huyghens' 
principle) in the present model. So, in a sense, the model also describes the return 
to (local) equilibrium. These results show that the arrow of time manifests itself 
already in the macroscopic space-time patterns of local equilibrium states. 

Of particular interest is the situation where ui describes a non-trivial macro- 
scopic hydrodynamic flow in some timelike cone, i.e. the local rest systems of u, 
which can be determined by the admissible macro-observables f3 \— > M e (f3) = e/3, 
e being positive lightlike, vary from point to point. In this case the arguments 
given in the preceding proposition imply that the function x h- > ui(M e )(x) has to 
be of the form 

u{M e )(x) = c w e(x - x w ), (5.3) 

where c w is some non-zero (say, positive) constant and x w some constant vector. If 
an observer finds such a behaviour in some bounded space-time region, he would 
infer that the maximal timelike cone where u can be locally in equilibrium is the 
lightcone V + + x u since u(M e )(x) has to be non-negative. Anticipating that uj 
has this property, he would also conclude that the functions x \— > u(M e )(x) vanish 
at the apex x^ of this cone, hence the temperature of the state tends to infinity 
there. Moreover, if x approaches any point on the ray x w + R + e at the boundary 
of the cone, the time axes of the respective local rest systems of the state bend 
towards the direction of e, indicating that there is a dominant flow of particles 
originating from x u . This can be made more explicit by looking at the particle 
density x i— > u(N p )(x) which can be shown to diverge if x — x^ and p become 
parallel. So the observer would be led to interpret the underlying state as the 
result of a hot bang which has occurred at x u and has quickly approached local 
equilibrium thereafter. Of course, all his conclusions depend on the hypothesis 
that the partial state which he observes is locally in equilibrium in all of V + + x^. 
Whether this is really the case can evidently not be decided by observations made 
in any finite space-time region. 

6 Summary and outlook 

In the present investigation we have applied the general formalism for the charac- 
terization and analysis of local equilibrium states, established in [T], to a simple 
model. There are two elemental ideas underlying this approach: (a) local equi- 
librium states cannot be distinguished from global ones by observations made in 
small space-time regions and (b) local measurements admit a macroscopic inter- 
pretation in all global equilibrium states. The combination of these two basic 
facts made it possible to interpret within the microscopic setting the space-time 
patterns of local equilibrium states in macroscopic terms. 

The results of our investigation show that this approach not only is physically 
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meaningful, but it also provides a convenient framework for concrete computations. 
In the case at hand, we have been able to establish equations for the space- 
time evolution of macroscopic properties of local equilibrium states, such as the 
particle density, without having to enter into the difficult dynamical question of 
how a system is driven towards local equilibrium. No a priori hypotheses, such 
as Boltzmann's assumption of molecular chaos, were necessary to establish these 
results. Instead, the intrinsic characterization (a) of local equilibrium was sufficient 
for their derivation. 

The present approach revealed also interesting constraints on the possible macro- 
scopic space-time structure of local equilibrium states. It turned out that local 
equilibrium can exist only in certain specific space-time regions which generically 
have the form of lightcones. Such states can therefore be used to distinguish be- 
tween future and past. Moreover, typically there appear singularities at the apices 
of these cones which have the interpretation of a hot bang. To the best of our 
knowledge, these results provide the first example of a "singularity theorem" in a 
microscopic setting. 

The present model of a massless free field exhibits a greatly simplifying feature, 
however: it is scale invariant, i.e. physics looks the same at microscopic and macro- 
scopic scales. For a study of the effects of an inherent scale, one may consider the 
example of massive free field theory. Indeed, there appear some differences. It is 
not difficult to show that, in the massive case, non-trivial local equilibrium states 
having locally a sharp temperature do not exist, their temperature support is typ- 
ically all of V+. This suggests to enlarge the space C of thermal reference states 
accordingly. In this extended framework it is still possible to establish evolution 
equations for the thermal equilibrium states, in complete analogy to the present 
results, and to analyze their space-time behaviour. These results will be put on 
record elsewhere. 

It seems attractive to apply, in a further step, the general formalism in pQ to 
free fields on arbitrary globally hyperbolic space-times. In particular, the shape of 
regions admitting local equilibrium is of interest there as well. For maximally sym- 
metric spaces there should appear no new conceptual difficulties in such an analysis 
since the basic ingredients, local observables and "global" equilibrium states for 
geodesic observers (i.e. states which are in equilibrium in the causal closure of the 
respective geodesies), do exist in these cases. But in generic spacetimes the situ- 
ation is less clear in view of the fact that such global equilibrium states need not 
exist. Thus what is required there is a substitute for these states which allows one 
to identify the local equilibrium situations. For the solution of this problem the 
generally covariant framework of quantum field theories, proposed in [S] , seems to 
be of relevance since it allows the unified treatment of a given theory on arbitrary 
spacetimes. So one may be able to use the existence of global equilibrium states on 
one spacetime to interpret the local equilibrium properties of states on the others. 

We mention as an aside that the present approach to the characterization of 
equilibrium states may be regarded as a refinement of the ideas underlying the 
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Hadamard condition jHj which is of fundamental importance in this setting. Yet 
whereas the Hadamard condition has been invented to characterize locally the 
folia of all physical states, the condition of local equilibrium is designed to identify 
within these folia states having a specific physical interpretation. Thus the present 
ideas seem to be of relevance also for the interpretation of states on arbitrary 
space-time manifolds. 

Returning to the case of Minkowski space, a central issue is the treatment of 
interacting models where one may study the intriguing question of the approach 
to equilibrium. In order to indicate the problems appearing in this context, let us 
consider first a state uo which is locally thermal on a space of observables generated 
by balanced derivatives 8^(x) of the normal ordered squares of the underlying 
fields [Ullj, similarly to the observables ()3.4|) used in the present investigation. In 
the interacting case, these observables no longer satisfy the equations ()4.3|) and 
(J4.4)) . which were the basis for the present analysis, there appear additional "source 
terms" a^(x), t^(x), 

d v e v ^{x) = a^(x), ne^(x) = t»(x). (6.1) 

Denoting by O' 1 the macro-observables corresponding to 8^(x), one thus obtains 
equations for their evolution in state u of the form 

d v u(e v f i )(x)=u(a' t ix)), Dui&^ix) =u(tI*(x)). (6.2) 

The terms appearing on the right hand side of these equations resemble the collision 
terms in transport equations 8\. In particular, they vanish if uj is sufficiently close 
to equilibrium at x since a^^x), r^(x) have zero expectation values in all thermal 
reference states due to their space-time invariance. 

In order to establish the approach to (local) equilibrium one thus has to show 
that these collision terms vanish asymptotically. The present results suggest that, 
for a given state u, one may expect such a behaviour at best in some lightcone 
and, because of the PCT theorem, this cone may be either future or past directed. 
So one has to identify the convex subsets of states exhibiting the desired behaviour 
for given time direction. Mixtures of the two types of states may in general not 
exhibit an ergodic behaviour. 

Next, there is the problem of local thermalization of a given state uo on the space 
of observables 8^(x), which was taken for granted in the preceding discussion. Its 
solution requires a proof that there are asymptotic bounds of the form pQ 

where || • \\b denotes any one of the seminorms on the macro-observables, intro- 
duced in Sec. 2. Choosing B sufficiently large, it seems feasible to establish such a 
result without major difficulties for finite dimensional spaces of observables 6 fl (x) 
whose corresponding thermal functions 0^ are linearly independent. Yet if there 
exist relations between these functions of the form Y2 c v = which do not have 
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a counterpart at the microscopic level, one has to understand the mechanism by 
which the expectation values of the corresponding combinations X] c a* of lo- 

cal observables vanish asymptotically in the state to. It is of interest in this context 
that there exists a large family of states which may be regarded as perturbations 
of global equilibrium states and which satisfy the above inequality pQ. 

Once local equilibrium has been reached in a state u, either exactly or in the 
sense that the collision terms in ()6.2j) may effectively be neglected, one arrives at 
the same conservation laws and evolution equations for the macro-observables 
in the state u as in free field theory, 

d v u(Q v,i )(x)=0 1 Dcuie^ix) = 0. (6.4) 

These equations thus have a universal character. Yet in order to be able to inter- 
pret them in physical terms it is necessary to determine the underlying thermal 
functions (3 i— > /x (/3), which are model dependent. This requires a computation 
of the thermal two-point functions of the underlying fields. The general form of 
these functions is known p] and there exist also various methods for their pertur- 
bative computation in models [TTHlllj. cf. also the novel approach in ^2]- It seems 
therefore possible to determine the macroscopic space-time patterns of local equi- 
librium states in the presence of interaction by similar methods as in the present 
investigation. We hope to return to this interesting issue elsewhere. 

Let us finally mention that the basic ideas underlying the approach in pQ can 
also be applied to non-relativistic theories, such as spin systems. Of course, the 
concept of local observable, as used in the present analysis, is no longer meaningful 
there, one has to replace it by suitable subspaces of the algebra of all observables. 
A general characterization of these subspaces is still an open problem. But a 
preliminary study ^3] of the XX-model indicates that certain macroscopic aspects 
of local equilibrium states, such as the existence of evolution equations, can still 
be established in such a setting, at least in simple cases. So there is evidence that 
the general approach to the analysis of local equilibrium states, proposed in [I], 
provides an efficient framework for the analysis of this important issue. 
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